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We report on a novel approach to cavity-ring-down spectroscopy with the sample gas in saturated- 
absorption regime. This technique allows to decouple and simultaneously retrieve empty-cavity 
background and absorption signal, by means of a theoretical model that we developed and tested. 
The high sensitivity and frequency precision for spectroscopic applications are exploited to measure, 
for the first time, the hyperfine structure of an excited vibrational state of '^^0^^C^''0 in natural 
abundance with an accuracy of a few parts in 10~^^. 
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The availability of a molecular-spectroscopy technique, 
able to combine the ultimate performance in terms of sen- 
sitivity, resolution and frequency accuracy, can be cru- 
cial in many fundamental physical measurements. In- 
deed, profiting from the strength and ease of satura- 
tion of many mid-IR ro- vibrational transitions, this tech- 
nique could provide new insights in elusive quantum- 
mechanical effects encoded in molecules, such as: parity 
violation due to weak interactions [T]-[5], violation of the 
Fermi/Bose statistics and/or the symmetrization postu- 
late of quantum mechanics [IHB], time variation of fun- 
damental physical constants, e.g. the proton-to-electron 
mass ratio[71|5j. Such a technique could also represent a 
major step forward in trace-gas sensing. 

Concerning sensitivity, cavity ring-down (CRD) 
spectroscopy jOj [lOj has proven to be a good candidate 
technique to directly provide a quantitative measurement 
of gas absorption coefficient with a simple experimen- 
tal set-up. In principle, it is not limited by amplitude 
noise of the laser source, but only by detection shot 
noise. However, variations of the empty-cavity decay 
rate always prevent to achieve this ultimate limit and 
to average measurements over long times. The empty- 
cavity background could be subtracted by quickly switch- 
ing the radiation frequency between nearby longitudinal 
cavity modes. Nevertheless, this method cannot be ef- 
fective, since we observed that decay times of different 
cavity modes are affected by uncorrelated fluctuations. 
Other techniques (e.g. CRD heterodyne spectroscopy[TT] 
and NICE-OHMS p2]) are more sensitive than standard 
CRD, but they all are more complex and require fast and 
sensitive detectors, generally unavailable in the mid IR. 

In this work we propose a new spectroscopic technique, 
namely saturated-absorption cavity ring-down (SCaR), 
that improves the CRD sensitivity. We show that the 
decrease of the saturation level during each SCaR event 
makes our technique very effective in identifying and de- 
coupling any variation of the empty-cavity decay rate. 
Saturation effects have already been observed in CRD 
experiments and, in some conditions, they induce a de- 



viation of the ring-down signal from the perfectly expo- 
nential behaviour. Whenever non-linear absorption oc- 
curred in CRD spectroscopy, either it was considered as 
a disadvantage to be avoided [TSl [H], or it was only used 
to get Lamb dips as frequency markers for line centers, 
without a detailed treatment of saturation effects[15]. In 
some continuous-wave CRD experiments non-linear ef- 
fects were observed and discussed in two specific regimes. 
When the ring-down time is much shorter than the pop- 
ulation relaxation time, absorption loss keeps constant 
at the saturated value existing at the CRD start time 
and the decay will be simply exponential|16j . In the op- 
posite regime (the so-called "adiabatic" approximation), 
the saturation level of the optical transitions reaches the 
steady state at each point of the ring-down profile and 
a more complete and complex theoretical model would 
be needed to fit the non-exponential decay curves|17j. 
Theoretical mean-field analyses on dynamic absorption 
saturation in pulsed CRD were performed a few years 
ago, either considering inhomogeneous broadening [TB] or 
not[19i. 

We developed and tested a new model which is very 
effective in exploiting the SCaR spectroscopic technique. 
We are assuming that the gas interacts with intra-cavity 
radiation in a TEMqo mode with a time-dependent inten- 
sity / and power P given by the following expressions: 

/(p,^) = /o(^)e-2(^)^ P{t) = ^-^m (1) 

where p = \/ x'^ + y'^ is the radial coordinate, /o(i) = 
I{p — 0, t) is the peak intensity on the cavity axis z and 
w is the beam waist, assumed to be constant along z. We 
also define several physical quantities: 
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where A is the Einstein coefficient of the transition, k is 
the wave vector and F is the full width at half maximum 
(FWHM) of the transition. We are neglecting the effects 
of the standing-wave light field inside the cavity and we 
are spatially averaging the different saturation levels of 
molecules interacting with light in node and anti-node 
positions. In the presence of inhomogeneous broaden- 
ing due to a thermal Gaussian distribution of molecular 
velocities, the absorption coefficient is affected by satu- 
ration and obeys the following equation [5D]: 



^l+I{p,t)/Is 
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where is the non-saturated value of a. Power atten- 
uation due to gas absorption, along the z axis, can be 
expressed as: 
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where the spatial integration over the transverse beam 
profile corresponds to a "local" approximation. That is 
correct at pressure regimes where the diffusion time of 
molecules is longer than the population relaxation time 
of molecular levels. Let us combine this continuous loss 
mechanism (following the Beer-Lambert law) with the 
discrete mirror losses (2 each cavity round trip), and de- 
fine the decay rates 7c = c(l — R)/l and 7g — cao? where 
R is the mirror reflectivity. Then, by using the speed-of- 
light relation ^ = we get the following rate equa- 
tion for the saturation parameter: 



dG 

~dt 



(t) = -7cG(t)-7, 
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G(0) = Go (5) 



where Gq = Pq/Ps is the saturation parameter value 
when the SCaR event starts, triggered by a threshold 
set on the signal detected in transmission. An analyt- 
ical solution of this equation can be found in the form 
t = t{G) that, unfortunately, cannot be inverted explic- 
itly in the form G — G{t). Instead, we chose to perform 
a numerical integration of Eqs. ([S]) with the 4*''-order 
Runge-Kutta algorithm (RK4) within the fitting proce- 
dure itself. Since 7g < 7c and the dynamic range of the 
decay curve exceeds 4 decades, in order to get a better 
computing precision let us factorize G as follows: 



G(i) = Goe-^=*/(i;Go,7c,7s) 



(6) 



With a direct substitution in Eq. ([s]), we can easily find 
that the function / obeys the following differential equa- 
tion: 
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The experimental set-up [H], is based on a difference- 
frequency-generated CW coherent source widely tunable 
in the mid IR, with the near-IR pump/signal lasers 
phase-locked one another through a fs Ti:sapphire op- 
tical frequency comb (OFC)[32]. The 1-m-long cavity 
is formed by 2 high-reflectivity mirrors with 6-m ra- 
dius of curvature and optical losses of 440 ppm around 
2340 cm~^. With this set-up we performed several spec- 
troscopic measurements to test both sensitivity and reso- 
lution using the newly developed model. Each spectrum 
is recorded by stepwise scanning the absolute frequency 
of the near-IR signal (and thus of the mid-IR idler). At 
each IR frequency the cavity length is dithered at about 
1 kHz rate across the resonance condition. When a fixed 
threshold level is reached during each cavity build-up, 
the IR frequency is rapidly switched off-resonance. Sev- 
eral CRD events are detected by a liquid-N2-cooled InSb 
photo-diode aligned with the cavity-transmitted light, ac- 
quired by a 18-bit digitizing oscilloscope with a sampling 
rate of 10 MS/s, averaged by a Lab VIEW acquisition 
program and processed on-line by a FORTRAN fitting 
routine. 

In Fig. [l] SCaR measurements are compared with the 
model we have developed. The experimental data in this 
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FIG. 1: Comparison between experimental data and theoret- 
ical model. Experimental data points for the case Go = 50, 
the fit curve (in red) and the residuals are plotted with left 
scale (in black). Saturated decay functions /, simulated for 
different values of Go (labeling curves), are plotted with right 
scale (in blue). 

figure are the average of 3072 decay signals measured at 
a fixed frequency near the absorption peak of the tran- 
sition of Fig. [2] over a 4-s time interval. The discrete 
noise values in the tail of the SCaR signal are due to 
the resolution of the digitizing oscilloscope. Almost flat 
residuals witness the validity of the theoretical model. A 
rough and intuitive explanation for the intrinsic ability of 
this technique to distinguish between empty-cavity and 
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gas-induced decay rates can be given as follows. Three 
consecutive intervals can be recognized in the SCaR sig- 
nal: zones A, B and C carry information respectively on 
7c, 7c + Ig and the detector offset due to thermal back- 
ground (whose fitted value has been subtracted from the 
signal in Fig. [I]) . The A-B transition is marked by a slope 
change in the SCaR signal, which needs the condition 
Go ^ 1 to be well observable, as it is more evident in the 
/ curves. As a consequence, a large detection dynamics is 
needed to measure the transition from strong-saturation 
to linear-absorption regime in the SCaR decay, before 
it falls below the noise level (B-C transition). The best 
choice for Gq should give such three zones with similar 
durations. 

Fig. [2] shows a low-pressure, Doppler-limited absorp- 
tion spectrum of a ro-vibrational transition of ^^C^^02 
recorded with SCaR. The free parameters in the fitting 
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FIG. 2: Doppler-limited absorption spectrum of the (03'^ 1 — 
03^0) R(50) transition of ^^C^^02 (line-center frequency v = 
2343.663 cm"^ line-strength 5* = 7.86 • lO"^"* cm), recorded 
with absorption length L = 1 m, temperature T ~ 296 K and 
pressure p = 50 (2) ^bar. The decay rate 7g (7c) is plotted 
with black (blue) scale at the left (right). The red curve fits 
7g experimental data points to a Gaussian (parameters in the 
inset and residuals at the bottom). 

routine were ^cii 7ci 7gi while was kept fixed, because it 
is frequency independent. The non-zero background ab- 
sorption is due to Lorentzian wings of nearby transitions. 
This spectrum was recorded in about 12 minutes by sep- 
arately averaging the 7g and 7c values resulting from the 
SCaR fitting routine. Since the IR frequency was scanned 
forward and backward in 10-MHz steps, no sub-Doppler 
features are visible at line center. The measured peak 
absorption rate H = 56.466 (33) ms^^ is consistent with 
the value 63 ms~^ calculated from the line-strength S 
reported by the HITRAN database|23l |24|, within the 
5-10% uncertainty on S and the 4% uncertainty on mea- 
sured pressure. A minimum detectable absorption coef- 
ficient amin = 1.1 ■ 10~^ cm~^ can be derived dividing 



by c the uncertainty value of H yielded by the fitting 
routine. When comparing Gaussian fit residuals of 7g 
in Fig. [2] with the measured 7c, it is evident that our 
technique effectively makes the measured gas absorption 
almost unaffected by any variations of the empty-cavity 
decay rate. Indeed, the fluctuations of the former curve 
are about a factor of 20 lower than the latter one and 
only a slight residual correlation is observed. This fac- 
tor also quantifies the sensitivity improvement over stan- 
dard CRD, that can only measure an overall decay rate 
7 = 7g -t- 7c and is inevitably limited by any background 
fluctuation contained in 7c. As mentioned above, this 
key point would enable a further gain of sensitivity with 
measurements averaged over longer times. 

To test the resolution, we performed sub-Doppler mea- 
surements on low-pressure ^^O^^C^^O at natural abun- 
dance (7.5 • 10"'*). Fig. [3] shows the resolved hyper- 
fine structure of the (00" 1 - 00"0) R(0) transition of 
^^O^^C^^O, due to interaction between the ^^O electric 
quadrupole (which is non-null for a nuclear spin 7 = 5/2) 
and the electric field gradient at the nucleus position[5S]. 
This spectrum was recorded in about 3 hours with 11 
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FIG. 3: Sub-Doppler spectrum of the (00^1 -00"0) R(0) tran- 
sition of ^^O^^C^'^O {u = 2340.765 cm-\ S = 1.25- 10"^^ cm), 
recorded with L = 1 m, T ~ 296 K and p ~ 2 ^bar. The 
dimensionless quantity D displaying the Lamb-dip features 
is defined as D = GoPa/Po- The experimental data and a 
multi-Lorentzian curve fitting the 3 Lamb dips and the cor- 
responding 3 cross-overs, are plotted. 

forward/backward frequency scans in 20-kHz steps. The 
experimental conditions, in this case, did not satisfy the 
local assumption of our model, since at such a low pres- 
sure the mean free path is larger than the beam diam- 
eter and the population relaxation time is dominated 
by the transit time tj. Therefore we dropped the spa- 
tial integration over the transverse beam profile of Eq. |4] 
and made a mean-field approximation. Moreover, the 
adiabatic approximation Tt <^ l/7c is at the validity 
edge. Nevertheless, this simplified model resolved well 
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the hyperfine structure. In this case the free parame- 
ters were Po,jc,D/Ps = Gq/Po while jg was kept fixed, 
because the gas absorption is almost constant at the 
Doppler peak value during the narrow frequency scan. 
We report in Table |T] the fit results for the line centers. 
Systematic errors due to the 10^^^ relative accuracy of 
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70174 358 594.9 (6.2) 


0.207(11) 





70174 357409.8 (3.0) 


0.345 (9) 


-fl 


70174 358 229.8 (3.6) 


0.448 (13) 



TABLE I: Measured absolute frequencies uaf of the hyperfine 
triplet with 1-a uncertainties and relative intensities. 

the absolute frequency scale (70 Hz) and to pressure- 
shift effects (-200 Hz) are negligible in the error budget, 
when compared with the statistical uncertainty (a few 
kHz) of the fitting routine. The relative intensities of the 
triplet components are nearly consistent with the theoret- 
ical 2:3:4 ratios yielded by the J, J, F quantum numbers 
involved|25). The fitted FWHM of each Lorentzian dip 
is r = 217.1 (5.5) kHz. The contributions to FWHM in 
our experimental conditions (w = 1.5 mm, p w 2 /ibar) 
are natural (22 Hz), pressure (14 kHz), and interaction- 
time broadening. The latter one is due both to transit 
time of molecules through the beam (55 kHz) and to 
mean lifetime of photons during a CRD event (21 kHz). 
We believe that the discrepancy with the larger mea- 
sured value of r can be ascribed to power broadening 
effects, due to the high average value of G during the 
SCaR process. The line-center frequency Vc and eqQ pa- 
rameters are obtained by fitting the hyperfine triplet to 
the theoretical expression depending on the /, J, F and 
I' , J' , F' quantum numbers. The measured line-center 
frequency i^c — 70 174 358 037.3 (3.9) kHz is consistent 
with the value 70174 368 MHz reported by HITRAN, 
within the declared 3-30 MHz uncertainty. Thus, we 
improved the frequency accuracy by more than 3 orders 
of magnitude. The measured electric-quadrupole cou- 
pling constant in the ro- vibrational state |00°1; 5/2, 1, F) 
is eqQ ~ —3.927 (30) MHz. A much earlier microwave 
measurement pS] of eqQ in the ground vibrational state 
was performed on a ^^0^'^C^^O isotopically enriched sam- 
ple, with an uncertainty about a factor of 3 worse than 
our measurement in a excited vibrational state. 

In conclusion, we have shown that SCaR spectroscopy 
is an effective novel approach to CRD. In particular, the 
saturation regime has proven to be generally useful for 
background reduction, whereas it also provides detection 
of sub-Doppler spectroscopic features, when higher reso- 
lution is needed. Finally, we have shown that frequency- 
comb-assisted SCaR spectroscopy can provide absolute 
frequency measurements of molecular transitions. 
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